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We investigate in a phenomenological way, direct CP violation in the hadronic 
^p Hi decays B^' — > tt + tt~ K^'° where the effect of p — uo mixing is included. If 

Oh, (the effective parameter associated with factorization) is constrained using the most 

recent experimental branching ratios (to p° K° , , p^ K° , p K^ and loK^) from 

the BABAR, BELLE and CLEO Collaborations, we get a maximum CP violating 
asymmetry, a max , in the range —25% to +49% for B~ — > it + it~K~ and —24% to 
+55% for B° — > 7r + 7r~ K°. We also find that CP violation is strongly dependent on 
the Cabibbo-Kobayashi-Maskawa matrix elements. Finally, we show that the sign 
of sin 5 is always positive in the allowed range of N% and hence, a measurement 
of direct CP violation in B^' Q — > tt^tt^ K ±,Q would remove the mod(7r) ambiguity 
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1 Introduction 



The study of CP violation in B decays is one of the most important aims for the B facto- 
ries. The relative large CP violating effects expected in B meson decays should provide 
efficient tests of the standard model through the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix. It is usually assumed that a nonzero imaginary phase angle, i], is responsible for 
the CP violating phenomena. This is why, in the past few years, numerous theoretical 
studies and experiments have been conducted in the B meson system [p], |2| in order to re- 
duce uncertainties in calculations (e.g. CKM matrix elements, hadronic matrix elements 
and nonfactorizable effects) and increase our understanding of CP violation within the 
standard model framework. 

Direct CP violating asymmetries in B decays occur through the interference of at least 
two amplitudes with different weak phase <fi and strong phase 5. In order to extract the 
weak phase (which is determined by the CKM matrix elements) through the measurement 
of a CP violating asymmetry, one must know the strong phase 5 and this is usually not 
well determined. In addition, in order to have a large signal, we have to appeal to some 
phenomenological mechanism to obtain a large 5. The charge symmetry violating mixing 
between p° and uj can be extremely important in this regard. In particular, it can lead to 
a large CP violation in B decays, such as B ±,Q — > p°(uj)K ±, ° — > ty+'k^K^ , because the 
strong phase passes through 90° at the uj resonance f|, [5| . 

We have collected the latest data for b to s transitions concentrating on the CLEO, 
BABAR and BELLE branching ratio results in our approach. The aim of the present 
work is multiple. The main one is to constrain the CP violating calculation in _B ±,CI — > 
p°(uj)K ±, ° — > 7i + 7i^K ±, °, including p — uj mixing and using the most recent experimental 
data for B — > pK decays. The second one is to extract consistent constraints for B decays 
into p{PS) where PS can be either tt or K. In order to extract the strong phase 5, we 
shall use the factorization approach, in which the hadronic matrix elements of operators 
are saturated by vacuum intermediate states. Moreover, we approximate non-factorizable 
effects by introducing an effective number of colors, . 

In this paper we investigate five phenomenological models with different weak form factors 
and determine the CP violating asymmetry, a, for B ±,Q — > p°(uj)K ±,(} — > -k + -k~ K^' in 
these models. We select models which are consistent with all the data and determine the 
allowed range for N^ ff (0.66(0.61) < A" c e// < 2.84(2.82)). Then, we study the sign of sin 5 
in this range of N%** for all these models. We also discuss the model dependence of our 
results in detail. 

The remainder of this paper is organized as it follows. In Section 2, we present the form 
of the effective Hamiltonian which is based on the operator product expansion, together 
with the values of the corresponding Wilson coefficients. In Section 3, we give the phe- 
nomenological formalism for the CP violating asymmetry in decay processes including 
p — uj mixing, where all aspects of the calculation of direct CP violation, the CKM ma- 
trix, p — uj mixing, factorization and form factors are discussed in detail. In Section 4, 
we list all the numerical inputs which are needed for calculating the asymmetry, a, in 
B ±,Q — ► p°(uj)K ±, ° — ► 7i + 7T^K ± ' . Section 5 is devoted to results and discussions for these 
decays. In Section 6, we calculate branching ratios for decays such as B ±,Q —>■ p ±,0 K ± '° 
and 5 ± — > ujK ± as well, and present numerical results over the range of iVj" which is 
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allowed by experimental data from the CLEO, BABAR, and BELLE Collaborations. In 
the last section, we summarize our results and determine the allowed range of N^** which 
is consistent with data for both pit and pK decays. Uncertainties in our approach and 
conclusions are also discussed in this section. 



2 The effective Hamiltonian 

2.1 Operator product expansion 

Operator product expansion (OPE) |J is a useful tool introduced to analyze the weak 
interaction of quarks. Defining the decay amplitude A[M — > F) as 

A(M - F) cc Ci(ji)(F\Oi(ji)\M) , (1) 



where Cj(/x) are the Wilson coefficients (see Section ^2|) and Oi(/x) the operators given by 



the OPE, one sees that OPE separates the calculation of the amplitude, A(M — > F), into 
two distinct physical regimes. One is related to hard or short- distance physics, represented 
by Ci(p) and calculated by a perturbative approach. The other is the soft or long-distance 
regime. This part must be treated by non-perturbative approaches such as the 1/N 
expansion [0], QCD sum rules || or hadronic sum rules. 

The operators, Oj, are local operators which can be written in the general form: 

O n = {qiT n iqj){qk^n2qi) , (2) 

where r n i and T n 2 denote a combination of gamma matrices and q the quark flavor. They 
should respect the Dirac structure, the color structure and the types of quarks relevant 
for the decay being studied. They can be divided into two classes according to topology: 
tree operators (Oi, 2 ), and penguin operators (0 3 to O w ). For tree contributions (W ± is 
exchanged), the Feynman diagram is shown Fig. [I]. The current-current operators related 
to the tree diagram are the following: 

0{ = q a J»(l - 75)^^7^(1 - 7 5 )& a , 

O s 2 = q la (l- 75)^(1 -75)6, (3) 

where a and (3 are the color indices. The penguin terms can be divided into two sets. The 
first is from the QCD penguin diagrams (gluons are exchanged) and the second is from 
the electroweak penguin diagrams (7 and Z° exchanged). The Feynman diagram for the 
QCD penguin diagram is shown in Fig. |2] and the corresponding operators are written as 
follows: 

03 = ^(1-75)^^(1-75)^, 

qi 

4 = q a l»(l - 75)&/3 i'fO- ~ 75)9^ » ( 4 ) 

q' 
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5 = q lfl (l - l5 )b <?V(1 + IsW , 

q' 

6 = q a ^(l - 75)^ J2 fa^ 1 + 75)5^ 



(5) 



where q' = u, d, s, c. Finally, the electroweak penguin operators arise from the two Feyn- 
man diagrams represented in Fig. |^ (Z, 7 exchanged from a quark line) and Fig. [| (Z, 7 
exchanged from the W line). They have the following expressions: 



7 = 


3 - n 


-75)^e g ,gV(l + 75)g' , 

q' 


o 8 = 




-75)&/3^ e 9'^ M ( 1+ ^ 5 )^ 


o 9 = 




-7 5 )&^vgY(i-75)g', 

q' 


010 = 




-7 5 )&/3^ e 9'^7 M (i -7 5 )g« 
9' 



(6) 



where e q > denotes the electric charge of q'. 
2.2 Wilson coefficients 

As we mentioned in the preceding subsection, the Wilson coefficients 0, Ci(ji), represent 
the physical contributions from scales higher than \i (the OPE describes physics for scales 
lower than Since QCD has the property of asymptotic freedom, they can be calculated 
in perturbation theory. The Wilson coefficients include contributions of all heavy particles, 
such as the top quark, the W bosons, and the charged Higgs boson. Usually, the scale \x 
is chosen to be of order 0{m h ) for B decays. Wilson coefficients have been calculated to 
the next-to-leading order (NLO). The evolution of C(/i) (the matrix that includes Ci(fj)) 
is given by, 

C{ii) = U{^M w )C{M w ) , (7) 
where U(/i, M w ) is the QCD evolution matrix: 



Ufa, M w ) 



47T 



a s {M ] 



Wi 



47T 



J 



with J the matrix summarizing the next-to- leading order corrections and U°(fi, Mw) the 
evolution matrix in the leading-logarithm approximation. Since the strong interaction 
is independent of quark flavor, the C(fj) are the same for all B decays. At the scale 
fi = mb = 5 GeV, C(fi) take the values summarized in Table |l] 0, [TTH . 

To be consistent, the matrix elements of the operators, Oi, should also be renormalized 
to the one-loop order. This results in the effective Wilson coefficients, which satisfy 
the constraint, 



C i (m b )(O i {m b )) = C' i (O i ) 



tree 



(9) 
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where (Oi) ree are the matrix elements at the tree level. These matrix elements will be 
evaluated in the factorization approach. From Eq. (^), the relations between C[ and Cj 
are 



10, 11 



where, 



c[ = 


Ci , 


C> 


= c 2 , 


c 3 = 


C 3 - P s /3 , 




= C A + P S 


c 5 = 


C 5 - P s /3 , 




= C 6 + P S 


c> = 


C 7 + P ei 




— Cg , 


c 9 = 


Cg + P e , 


Cio 


= Cio > 


Ps 


= (a s /87r)C 2 (10/9 + G(m c , 







p p 



(a em /97r)(3C 1 + C 2 )(10/9 + G(m c ,/x, q 2 )) 



(10) 



and 



G(m c , fi, q 2 ) = 4 / dxx(x — l)ln 



x(l — x)q 2 



(12) 



Here q 2 is the typical momentum transfer of the gluon or photon in the penguin diagrams 
and G(m c , fi,q 2 ) has the following explicit expression [|T^1 , 



3?e G 



3m G 




1-4^ 



4^fln— 
9 l-Jl-4 m 



(13) 



Based on simple arguments at the quark level, the value of q 2 is chosen in the range 
0.3 < q 2 /m 2 < 0.5 || ||. From Eqs. (p~Q| - p~3]) we can obtain numerical values for C[. 
These values are listed in Table where we have taken a s {m z ) = 0.112, a em {m h ) = 
1/132.2, m b = 5 GeV, and m c = 1.35 GeV. 



2.3 Effective Hamiltonian 

In any phenomenological treatment of the weak decays of hadrons, the starting point is 
the weak effective Hamiltonian at low energy Ul3| . It is obtained by integrating out the 
heavy fields (e.g. the top quark, W and Z bosons) from the standard model Lagrangian. 
It can be written as, ^ 

H eff = -| Yl VckmC^O^) , (14) 



where Gp is the Fermi constant, Vckm is the CKM matrix element (see Section |3~TD , C;(/z) 



are the Wilson coefficients (see Section pT2|) , 0i(/x) are the operators from the operator 



product expansion (see Section |2.1| ), and fx represents the renormalization scale. We 



4 



emphasize that the amplitude corresponding to the effective Hamiltonian for a given decay 
is independent of the scale \l. In the present case, since we analyze direct CP violation 
in B decays, we take into account both tree and penguin diagrams. For the penguin 
diagrams, we include all operators O3 to 0\q. Therefore, the effective Hamiltonian used 
will be, 

,ab=i Gf 



n 



eff 



MB PV) = «| 



V2 



10 



+ H.C. 



VutV^dOl + C 2 O s 2 ) - V tb V* CiO 

i=3 

and consequently, the decay amplitude can be expressed as follows, 

V ub V: s (Cx{PV\0{\B) + C 2 {PV\O s 2 \B))- 

10 

VaV^C^PVmB) 



(15) 



i=3 



+ H.c. , (16) 



where (PV\Oi\B) are the hadronic matrix elements. They describe the transition between 
the initial state and the final state for scales lower than /1 and include, up to now, the 
main uncertainties in the calculation since they involve non-perturbative effects. 



CP violation in B 



±,o 



Direct CP violation in a decay process requires that the two CP conjugate decay processes 
have different absolute values for their amplitudes |TJ]. Let us start from the usual 
definition of asymmetry, 



which gives 



a(B -> F) 



a{B -> F) 



T(B -»• F) - T{B -> F) 
~ T(B -* F) + T(B -»• F) ' 

\A(B ^ F)\ 2 - \A(B ^ F)\ 2 



(17) 



\A(B -> F)| 2 + \A(B -> F)| 2 ' v ' 

where A [B — > F) is the amplitude for the considered decay, which in general can be 
written as A(B F) = |v4i|e i,5l+i<?il + \A 2 \e i&2+i<l>2 . Hence one gets 



a{B -> F) 



— 2|v4i 


||A 


2 


sin(0i - (j) 2 ) sin(5i - 8 2 ) 




2 + 2|A!|| 


^2 


cos(0i - 2 ) cos(5i - 5 2 ) + 


A 2 


2 



(19) 



Therefore, in order to obtain direct CP violation, the CP asymmetry parameter a needs 
a strong phase difference, S\ — 5 2 , coming from the hadronic matrix and a weak phase 
difference, 0i — 2 , coming from the CKM matrix. 



3.1 CKM matrix 

In phenomenological applications, the widely used CKM matrix parametrization is the 
Wolfenstein parametrization . In this approach, the four independent parameters are 
\,A,p and n. Then, by expanding each element of the matrix as a power series of the 
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parameter A = sin 9 C = 0.2209 (0 C is the Gell-Mann-Levy-Cabibbo angle), one gets (0(A 4 ) 
is neglected) 



/ 1 - |A 2 A AX 3 {p-i V )\ 

Vckm =\ -A 1 - \X 2 AX 2 , (20) 

\ AX 3 (1 -p-irf) -AX 2 1 / 

where rj plays the role of the CP- violating phase. In this parametrization, even though 
it is an approximation in A, the CKM matrix satisfies unitarity exactly, which means, 

Vckm ' Vckm = I = Vckm • V CKM . (21) 



3.2 p — uj mixing 



In the vector meson dominance model [16|, the photon propagator is dressed by coupling 



to vector mesons. From this, the p — uj mixing mechanism [17fl was developed. Let A be 
the amplitude for the decay B — > p°(u)K — > 7t + 7t~K, then one has, 



A = (Kir~7r + \H T \B) + (Kir-ir + \H p \B) 



[22] 



with H T and H p being the Hamiltonians for the tree and penguin operators. We can 
define the relative magnitude and phases between these two contributions as follows, 



A = (Kn-n + \H T \B)[l + re i& e i4> ] , 
A = (Kn + n-\H T \B)[l + re w e~^] 



(23) 



where 5 and <fi are strong and weak phases, respectively. The phase arises from the 
appropriate combination of CKM matrix elements, and = a.rg[(V t bV t *) / (V u bV* s )]. As a 
result, sin0 is equal to sin7 with 7 defined in the standard way [|18|| . The parameter, r, 
is the absolute value of the ratio of tree and penguin amplitudes: 



r = 



(p°(u)K\H p \B) 



(p°(uj)K\H T \B} 



(24) 



In order to obtain a large signal for direct CP violation, we need some mechanism to 
make both sin S and r large. We stress that p — uj mixing has the dual advantages that 
the strong phase difference is large (passing through 90° at the uj resonance) and well 
known [|], |5|. With this mechanism (see Fig. |5|), to first order in isospin violation, we have 
the following results when the invariant mass of is near the uj resonance mass, 



(Kir-ir+lH 1 \B) 
(Kir-TT + \H P \B) 



t + —t 



9 P fr 



+ —P P 

Sp 



(25) 



Here ty (V = p or uj) is the tree amplitude and py the penguin amplitude for producing 
a vector meson, V, g p is the coupling for p° — > tt + ti^ , n pw is the effective p — uj mixing 
amplitude, and sy is from the inverse propagator of the vector meson V, 



Sy = s — ml + im v T 



v 



(26) 
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with y/s being the invariant mass of the 7r + 7r~ pair. We stress that the direct coupling 
uj — > 7r + 7r _ is effectively absorbed into Yi pu} []19| , leading to the explicit s dependence of 
Hpuj. Making the expansion U^s) = U puj (ml) + (s — m^)n^(m^), the p — uj mixing 
parameters were determined in the fit of Gardner and O'Connell [pOfl : 3fte ~n. pu) {ml) = 
-3500 ± 300 MeV 2 , 3m %^{ml) = -300 ± 300 MeV 2 , and ti' pu) (ml) = 0.03 ± 0.04. In 
practice, the effect of the derivative term is negligible. ^From Eqs. ^) one has 



Defining 

— = r'e l ^ +<t>) , - = ae iS « , ^ = (3e i& ? , (28) 

h p h p Flu 

where 5 ai 5p and 5 q are strong phases (absorptive part). Substituting Eq. ( |28| ) into 
Eq. fl2~7|), one finds: 

re" = r'e iS * ^ t^*"*? , (29) 

where 



ae iS ° = f, (3e i5 ' 3 = b + ci, r'e iS " =d+ei , (30) 
and using Eq. (|29|) , we obtain the following result when a/s ~ m w : 



re 



C + iD 



(s-ml + f3ie II p J 2 + (/Sm II pw + 
Here C and .D are defined as: 



(31) 



C= (s-ml + We U^ld 



9fte Il pw + 6(s - ml) - cm u r, 



— e 



Il pa; + bm^T^ + c(s - m u 



+ (/3m Il pa; + m w r w ) < e 



3f?e n pw + 6(s - ml) - cm u T u 



+ d 



3m U pw + bm u Y u + c(s - m u 



, (32) 
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and 



D= (s-ml + f^e tl p ^)\e 



5te Tl pu) + d(s - ml) - cm^Tu 



+ d 



Sm I\. pu) + bm^T^ + c(s - m„ 



(f^m fl pu) + m/ u ) <^ d 



5?e n pw + b{s - ml) - cm u T u 



Q777. Yi pu) + bm^T^ + c(s - m, 2 



• (33) 



ae t5a , j3e tSl3 , and rV 5,7 will be calculated later. In order to get the CP violating asymmetry, 
a, sin0 and cos</> are needed, where is determined by the CKM matrix elements. In 
the Wolfenstein parametrization [[Tj|, the weak phase comes from [Vt&V^/V^V^J and one 
has for the decay B — > p°(uj)K, 

-v 



sin . 



COS( 



Vp 2 + v 2 



rj- 



(34) 



The values used for p and r\ will be discussed in Section 4.1 



3.3 Factorization 

With the Hamiltonian given in Eq. (p~5| ) (see Section |2.3| ), we are ready to evaluate 
the matrix elements for B ±,Q —>■ p°(u)K ±, °. In the factorization approximation |2jfl , 
either p°(uj) or K ±,a is generated by one current which has the appropriate quantum 
numbers in the Hamiltonian. For these decay processes, two kinds of matrix element 
products are involved after factorization (i.e. omitting Dirac matrices and color la- 
bels): (p^u^iuu^iK^Ksb)^ ) and (^ ± '°|(gig 2 )|0)(p°(cj) 1(^)1^°), where q x and 
q2 could be u, s or d. We will calculate them in several phenomenological quark models. 
The matrix elements for B — > X and B —>■ X* (where X and X* denote pseudoscalar and 



vector mesons, respectively) can be decomposed as follows [22 



(X\J,\B) = {p B + Px - ml k f x k\ F 1 (k 2 ) + ml k f x k,F {k 2 ) , (35) 



and, 



2a 



(X*\J p \B) = ■ e^ pa e*»p p B p x *V(k 2 ) + i{ e*(m B + m x *)A x {k 

m B + m x * - ' I 

fk ' k (P B + Px*)»A 2 {k 2 ) - € ^2m x * ■ k p A 3 (k 2 ) 



m B + m x * k A 

.€*■ k 



+ t—2m x *-k ll A (k 2 ) , (36) 
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where is the weak current, defined as J M = (?7 M (1 — 75)6 with q = u,d,s and 
k = pb — Px(x*)- tfj, is the polarization vector of X*. F and Fi are the form factors 
related to the transition 0~ — > CT, while Aq, Ai, A2, A3 and V are the form factors that 
describe the transition 0~ — ► Finally, in order to cancel the poles at q 2 = 0, the form 
factors respect the conditions: 



F 1 (0)=F (0), A 3 (0) = A (0) 
and they also satisfy the following relations: 



A 3 (k 2 



A,(k 2 



tub — mx* 



A 2 {k 2 



(37) 



(38) 



2mx* 2mx+ 

An argument for factorization has been given by Bjorken the heavy quark decays 
are very energetic, so the quark-antiquark pair in a meson in a final state moves very fast 
away from the localized weak interaction. The hadronization of the quark-antiquark pair 
occurs far away from the remaining quarks. Then, the meson can be factorized out and 
the interaction between the quark pair in the meson and the remaining quark should be 
tiny. 

In the evaluation of matrix elements, the effective number of colors, Nf;** , enters 
through a Fierz transformation. In general, for operator Oj, one can write, 



1 



(N, 



g + 6 >with * = !,••• ,10 



(39) 



where & describes non-factorizable effects. We assume & is universal for all the operators 
Oi. We also ignore the final state interactions (FSI). After factorization, and using the 
decomposition in Eqs. (|3lS|, BBI), one obtains for the process B° — > p (u)K°, 



t P = m B \p p \(C[ + —C' 2 )f p F l 



m 



(40) 



where f p is the p decay constant (and to simplify the formulas we use N c for N%** in Eqs. 
(40)-(50)). In the same way, we find t w = t p , so that, 

6<< 



ae"" a = 1 . 

After calculating the penguin operator contributions, one has 



r ' e iS q 



Pu 



(C[ + ^C' 2 )f p F l (m^ p ) 



v tb v* 



(41) 



(42) 



and 



m B \p P 

Pu 



ni(C' 7 + W C>) + (C 9 + — Cio) ) f P FMD+ 



N r 



+ 



K 



(m b + m d )(m d + m s ) _ 



f K A (m 2 K )} , (43) 
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where fx is the K decay constant. In Eqs. ((42], [43]), p^ has the following form: 



+ i(K + ^Q + (Q + 7 L C ; ))/,F ] 



m 2 K f K A (m 



_{m b + m d ){m d + m s 
and the CKM amplitude entering the 6 — > s transition is, 



(44) 



l l 



l l 



A2 y^Trf A 2 |sin/?| 



(45) 



with /? defined as the unitarity triangle as usual. Similarly, by applying the same formal- 
ism, one gets for the decay B~ — > p°(u)K~ , 



t p = m B 



1 



(C[ + F Ci)/ p Pi(m;) + (C5 + —C[)f K A (m 2 K 



In the same way, we find t w = t p , therefore one has again, 



ae 



JSa 



1 . 



(46) 



(47) 



The ratio between penguin and tree operator contributions, which involves CKM matrix 
elements, is given by, 



r'e iS " 



Pu 



v, lh v* 



and finally, 



(48) 



0e 



m B \p 

Pu 



r]) {C' A + —C^f K A Q {ml) 



+ I (Vr + ±-C>) + (C> + ±-C[ )^ fpFtiml) + (C' w + ±-C>)f K A (ml) 

m 2 K f K A (m 2 K ) 



_{m u + m s ){m b + m u ) _ 



(49) 
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where the u penguin contribution, p w , is: 



P. = m B \p p \U((C' z + -Left + (C 5 + ^C',))f P Hrn 2 p ) 
+ \ ((C 7 + ^-C 8 ) + (C 9 + i-C%)) f p F x (m 2 p ) + ((C' A + i<r s ) + (C[ + l-CrSj f K A (m 2 K ) 



m 



K 



(m u + m s ){m b + m u ) 



f K A (m 2 K ) \ . (50) 



3.4 Form factors 

The form factors Fi(k 2 ) and Aj(k 2 ) depend on the inner structure of the hadrons. We 
will adopt here three different theoretical approaches. The first was proposed by Bauer, 
Stech, and Wirbel ]22| (BSW), who used the overlap integrals of wave functions in order to 
evaluate the meson- meson matrix elements of the corresponding current. The momentum 
dependence of the form factors is based on a single-pole ansatz. The second one was 
developed by Guo and Huang (GH) pi| . They modified the BSW model by using some 
wave functions described in the light-cone framework. The last model was given by 
Ball |25| and Ball and Braun |26|]. In this case, the form factors are calculated from QCD 
sum rules on the light-cone and leading twist contributions, radiative corrections, and 
5'L r (3)-breaking effects are included. Nevertheless, all these models use phenomenological 
form factors which are parametrized by making the nearest pole dominance assumption. 
The explicit k 2 dependence of the form factor is as |22|. 2l, p5l 



10: 



F 1 (k 2 



hi 



k 2 



or 



hi 



+ h ( 

m B \ 



k' 2 



2 • 



A (k 2 



A (k 2 



h 



Ao 



k 2 



h 



An 



k 2 



2 ' 



(51) 



where n = 1,2, rriA and mi are the pole masses associated with the transition current, 
hi and hA are the values of form factors at q 2 = 0, and di and foj (i — 0, 1) are parameters 
in the model of Ball. 



4 Numerical inputs 

4.1 CKM values 

In our numerical calculations we have several parameters: q 2 , N^*, and the CKM matrix 
elements in the Wolfenstein parametrization. As mentioned in Section 1233, the value of q 2 



is conventionally chosen to be in the range 0.3 < q 2 /m b 2 < 0.5. The CKM matrix, which 
should be determined from experimental data, is expressed in terms of the Wolfenstein 
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parameters, A, A, p, and r] [|15j. Here, we shall use the latest values [^8| which were 
extracted from charmless semileptonic B decays, (|V^&|), charm semileptonic B decays, 
(|Vc*|), s and d mass oscillations, Am s , Am^, and CP violation in the kaon system (ex), 
(p,rj). Hence, one has, 

A = 0.2237 , A = 0.8113 , 0.190 < p < 0.268 , 0.284 < r] < 0.366 . (52) 

These values respect the unitarity triangle as well (see also Table |3[). 

4.2 Quark masses 

The running quark masses are used in order to calculate the matrix elements of penguin 
operators. The quark mass is taken at the scale /i ~ m b in B decays. Therefore one 



has 29 



m u (p = m b ) = 2.3 MeV , rn d (p = m&) = 4.6 MeV , 

m s (n = m b ) = 90 MeV , m b (p = m b ) = 4.9 GeV , (53) 

which corresponds to m s (p = 1 GeV) = 140 MeV. For meson masses, we shall use the 
following values JT8] |: 



m B ± = 5.279 GeV , m B ° = 5.279 GeV , 

m K± = 0.493 GeV , m K o = 0.497 GeV , 

m^± = 0.139 GeV , m^o = 0.135 GeV , 

m p0 = 0.769 GeV , m u = 0.782 GeV . (54) 

4.3 Form factors and decay constants 

In Table (| we list the relevant form factor values at zero momentum transfer [ 22], p4"[ p5| , 



26] , pOl for the B — > K and B — > p transitions. The different models are defined as follows: 
models (1) and (3) are the BSW model where the q 2 dependence of the form factors is 
described by a single- and a double-pole ansatz, respectively. Models (2) and (4) are the 
GH model with the same momentum dependence as models (1) and (3). Finally, model 
(5) refers to the Ball model. We define the decay constants for pseudo-scalar (_/>) and 
vector (fv) mesons as usual by, 

CP(g)|?i7**7552|0) = ifpq^ , 
V2(V(q)\qx^\0) = f v m v e v , (55) 

with g M being the momentum of the pseudo-scalar meson, my and ey being the mass and 
polarization vector of the vector meson, respectively. Numerically, in our calculations, we 
take 



f K = 160 MeV , f p ~ f u = 221 MeV . (56) 

The p and us decay constants are very close and for simplification (without any conse- 
quences for results) we choose fp — f^. 
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5 Results and discussion 



We have investigated the CP violating asymmetry, a, for the two B decays: B° — > p°K° — > 
7r + 7r~^° and -B~ — > p°K~ — > 7r + 7r _ i^ _ . The results are shown in Figs. [] and [7| for B° — > 

1T + 1T-K , (a = [T(B° -> TT + TT-i? ) - -> 7T-7T + ir )]/[r(fi -> TT+TT' K°) + F(5° -> 

tt-tt+K )]), where k 2 /m 2 b = 0.3(0.5) and for N e Jf equal to 0.61,0.66,2.65,2.69,2.82 and 
2.84. Similarly, in Figs. |8| and [|, the CP violating asymmetry, a, (= \F(B~ — > 7t + tt~K~) — 
T(B + -> 7T-7r+ir + )]/[r( J B- -> tt+tt-^-) + T(B + -> 7^+^+ )] ) , is plotted for B~ -> 
7r + 7r~i^~, where k 2 /m 2 = 0.3(0.5) and for the same values of N^** previously applied for 
B° — > tt^tt^K . In our numerical calculations, we found that the CP violating parameter, 
a, reaches a maximum value, a max , when the invariant mass of the 7r + n~ is in the vicinity 
of the u resonance, for a fixed value of A^f ■ We have studied the model dependence of 
a with five models where different form factors have been applied. Numerical results for 
B° — > 7t + 7t~K° and B~ — > 7T + 7r _ i^ _ are listed in Tables ^| and|^, respectively. It appears 
that the form factor dependence of a for all models, and in both decays, is weaker than 
the dependence. 

For B° — > 7t + 7t~K°, we have determined the range of the maximum asymmetry param- 
eter, a max , when N^* varies between 0.66(0.61) and 2.84(2.82), in the case of k 2 /ml = 
0.3(0.5). The evaluation of a max gives allowed values from 37%(55%) to — 20%(— 24%) 
for the range of iV^" and CKM matrix elements indicated before. The sign of a max stays 
positive until iV^ reaches 2.7. If we look at the numerical results for the asymmetries 
(Table |5|), for N^J in = 0.66(0.61) and k 2 jm\ = 0.3(0.5), we find good agreement be- 
tween all the models, with a maximum asymmetry, a max , around 33%(45.6%) for the set 
(Pmax,Vmax), and around 26%(33.2%) for the set (p m in,Vmin)- The ratio between asym- 
metries associated with the upper and lower limits of (p,rj) is around 1.26(1.37). If we 
consider the maximum asymmetry parameter, a max , for N^J ax = 2.84(2.82), we observe a 
distinction between the models. Indeed, two classes of models appear: models (2) and (4) 
and models (1,3) and (5). For models (2) and (4), one has an asymmetry, a max , around 
— 6%(— 7%) and around — 9%(— 10%) for the upper and lower set of (p,r)), respectively. 
The ratio between them is around 1.50(1.42). For models (1,3) and (5), the maximum 
asymmetry is of order — 14.3%(— 16.3%) for (p m ax,Vmax) and around — 19.3%(— 23.0%) for 
(PmimVmin)- in this case, the ratio between asymmetries is around 1.34(1.41). 

The first reason why the maximum asymmetry, a max , can vary so much comes from 
the element V u b. The other CKM matrix elements Vtb, Vt s and V us , all proportional to A 
and A, are very well measured experimentally and thus do not interfere in our results. 
Only V u b, which contains the p and r\ parameters, provides large uncertainties, and thus, 
large variations for the maximum asymmetry. The second reason is the non-factorizable 
effects in the transition b — > s. It is well known that decays including a K meson (and 
therefore an s quark) carry more uncertainties than those involving only a tt meson (u,d 
quarks). If we look at the asymmetries at N^ in , all models give almost the same values, 
whereas at N%ff ax , we obtain different asymmetry values (with, moreover, a change of sign 
for the CP violating asymmetry). The CP asymmetry parameter is more sensitive to 
form factors at high values of N^* than at low values of N%** . It appears therefore that 
all of the models investigated can be divided in two classes, referring to the two classes of 
form factors. 
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For B~ — > 7T + 7r _ K _ , we have similarly investigated the CP violating asymmetry. 
The values of maximum asymmetry parameter, a max , for a range of N^f from 0.66(0.61) 
to 2.84(2.82), where k 2 jm\ = 0.3(0.5) and for the five models analyzed, are given in 
Table || We found that for this decay, the CP violating parameter, a, takes values 
around 49%(46%) to — 22%(— 25%) for the limiting CKM matrix values of p and rj defined 
before. Once again, the sign of the asymmetry parameter, a, is positive if the value of 
N^ff stays below 2.7. If we focus on N^f in equal to 0.66(0.61), models (1,2,3,4) and 
(5) give almost the same value which is around 46.6%(43.6%) for the maximum values 
of the CKM matrix elements. For the set {p m in,Vmin), the maximum asymmetry, a, is 
around 34.0%(33.8%). The ratio between asymmetry values taken at upper and lower 
limiting p and rj values is around 1.37(1.28). Let us have a look at the CP asymmetry 
values at N^ ax . As we observed for the decay B° — > n + 7r^K°, all models are separated 
into two distinct classes related to their form factors. For models (1, 3) and (5), the value 
of maximum asymmetry, a max , is around — 15.6%(— 17.6%) and around — 21%(— 23.6%) 
for the maximum and minimum values of set (p,rj), respectively. The calculated ratio is 
around 1.34(1.34), between these two asymmetries. As regards models (2) and (4), for 
the same set of (p,r]), one gets — 11.5%(— 13%) and — 17%(— 18%). In this case, one has 
1.47(1.38) for the ratio. The reasons for the differences between the maximum asymmetry 
parameter, a max , are the same as in the decay B° — > 7t + tt~K°. 

By analyzing the B decays, such as B° — > 7t + tt~K° and B~ — > tt + tt~K~ , we found that 
the CP violating asymmetry, a, depends on the CKM matrix elements, form factors and 
the effective parameter N^f (in order of increasing dependence). As regards the CKM 
matrix elements, the dependence through the element, V u t>, contributes to the asymmetry 
in the ratio between the u penguin contributions and the p tree contributions. It also 
appears that for the upper limit of set (p,r)), we get the higher value asymmetry, a, and 
vice versa. With regard to the form factors, the dependence at low values of N^** is very 
weak although the huge difference between the phenomenological form factors (models 
(2) and (4) and models (1, 3) and (5)) applied in our calculations. At high values of N*ff, 
the dependence becomes strong and then, the asymmetry appears very sensitive to form 
factors. For the effective parameter, N^*, (related to hadronic non-factorizable effects), 
our results show explicitly the dependence of the asymmetry parameter on it. Because 
of the energy carried by the quark s, intermediate states and final state interactions 
are not well taken into account and may explain this strong sensitivity. Finally, results 
obtained at k 2 jm\ = 0.3(0.5), also show renormalization effects of the Wilson coefficients 
involved in the weak effective hadronic Hamiltonian. For the ratio between asymmetries, 
results give an average value of order 1.36(1.40) for B° — > -k + 7t~K and 1.39(1.33) for 
B~ — > 7t + tt~K~. This ratio is mainly governed by the term 1/ sin/3, where the values of 
the angles a, (3 and 7 are listed in Table |3|. 

As a first conclusion on these numerical results, it is obvious that the dependence of 
the asymmetry on the effective parameter iVj" is dramatic and therefore it is absolutely 
necessary to more efficiently constrain its value, in order to use asymmetry, a, to determine 
the CKM parameters p and 77. We know that the effects of p — uj mixing only exist 
around uj resonance. Nevertheless, in Figs. |5], |7j, and |S], at small values of N^* , e.g. 
~ 0.6, the curves show large asymmetry values far away from u resonance, which is a 
priori unexpected. In fact, if we assume that nonfactorizable effects are not as important 
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as factorizable contributions, then N^* should be much bigger (see Eq. ([$9])). /,From 
previous analysis on some other B decays such as B — > Dir, B — > ujtt, and B — > ujK, it 
was found that N*" should be around 2 BT|. Therefore, although small values of N"x" 



are allowed by the experimental data we are considering in this paper, we expect that the 
value of iVg" cannot be so small with more accurate data. We have checked that when 
N^*' is larger than 1 the large CP asymmetries are confined in the uj resonance region. 
With a very small value of N^f*, nonfactorizable effects have been overestimated. This 
means that soft gluon exchanges between p°(ui) and K may affect p — uj mixing and hence 
lead to the large CP asymmetries in a region far away from uj resonance. However, when 
y/s is very far from uj resonance, the CP asymmetries go to zero as expected. 

In spite of the uncertainties discussed previously, the main effect of p — uj mixing in 
B — > 7t + tt~K is the removal of the ambiguity concerning the strong phase, sin S. In the 
b — > s transition, the weak phase in the rate asymmetry is proportional to sin 7 where 
7 = arg[— (V ts V£) / (V us V* b )]. Knowing the sign of sin 5, we are then able to determine 



the sign of sin 7 from a measurement of the asymmetry, a. In Figs. [T(] and [II], the 
value of sin 5 is plotted as a function of Nl*' for B° — ► n + Tf^K° and B~ — > n + Tf^K~, 
respectively. It appears, in both cases, when p — uj mixing mechanism is included, that 
the sign of sin 5 is positive, for all models studied, until Nf;** reaches 2.69(2.65) for both 
B~ — > it + tt~K~ and B° — > tt + tt~K° } when k 2 jm\ = 0.3(0.5). For values of N%H bigger 
than this limit, sin S becomes negative. At the same time, the sign of the asymmetry 
also changes. In Figs. [T2|b and [TB|b, the ratio of penguin to tree amplitudes is shown 



for B ±,Q — > tt + tt~ K^' , in the case of n pa; = (—3500,-300). The critical point around 
= 2.7, refers to the change of sign of sin 5. Clearly, we can use a measurement of 
the asymmetry, a, to eliminate the uncertainty mod(7r) which is usually involved in the 
determination of 7 (through sin 27). If we do not take into account p — uj mixing, the 
CP violating asymmetry, a, remains very small (just a few percent) in both decays. In 



Figs. and |TT] (for the evolution of sin S) and in Figs. |l2| a and |i3| a (for the evolution of 
penguin to tree amplitudes), for B ±,Q — > 7r + 7r we plot sin 5 and r when U pLU = (0, 0) 

-i.e. without p — uj mixing. There is a critical point at Nl'* = 1 (for B° — ► 7t + 7t~K ) 
and N^'J = 0.24 (for B~ — ► 7t + 7t~K~) for which the value of sin 5 is at its maximum and 
corresponds (for the same value of N^f), to the lowest value of r. The last results show 
the double effect of the p — uj mixing: the CP violating asymmetry increases and the sign 
of the strong phase S is determined. 

6 Branching ratios for B ±,Q — > p°K ±,Q 
6.1 Formalism 

With the factorized decay amplitudes, we can compute the decay rates by using the 
following expression 



T(B -> VP) - ' 



Truly 



A(B -> VP) 



(57) 
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where p p is the cm. momentum of the decay particles defined as, 



\Pp\ 



Vl m B ~ ( m i + m 2 ) 2 ][m B - (mi - m 2 ) 2 ] 



2m 



B 



(5f 



mi(m 2 ) is the mass of the vector (pseudo-scalar) V(P) particle, ey is the polarization 
vector and A(B — > VP) is the decay amplitude given by, 



A(B - VP) = % £ V^VmiB) 
v 2 i=l,l0 



(59) 



where the effective parameters, a*, which are involved in the decay amplitude, are the 
following combinations of effective Wilson coefficients: 



a 2j 



C' 2j + 



1 



n; 



rCL 



a 2j -i = C' 2i _ x + t^jjC'^, for j = 1 



n; 



5 . 



(60) 



All other variables in Eq. ( p^| ) have been introduced earlier. In the Quark Model, the 
diagram (Fig. [| top) gives the main contribution to the B — > decay. In our case, 
to be consistent, we should also take into account the p — uj mixing contribution (Fig. [| 
bottom) when we calculate the branching ratio, since we are working to the first order of 
isospin violation. The application is straightforward and we obtain the branching ratio 
for B -> p°K: 



BR(B -> p°K) 



G\\V P 



VjA T J ail a 2 )-V s p A p Ja^ 



a k nT B 

i V s T Al( ai ,a 2 )-V s p A p (a 3 , 



• O'lQ) 



n 



(sp - m 2 ) + im w Y u 



(61) 



In Eq. ( pip Crp is the Fermi constant, V B is the total width B decay, and a k is an integer 
related to the given decay. Ay and A p are the tree and penguin amplitudes which respect 
quark interactions in the B decay. V^' p (in Eq. (0)) or V" S T , V p (in Eq. (|BTP) represent 
the CKM matrix elements involved in the tree and penguin diagram, respectively: 



V 



\v ub v: 



for 



1,2 , and V p = WtbV+l for 



10 . 



(62) 



6.2 Calculational details 

In this section, we enumerate the theoretical decay amplitudes. We shall analyze five 
b into s transitions. Two of them involve p — uj mixing. These are B~ — > p°K~ and 
B° — > p°K°. Two other decays are B° — > p~K + and B~ — > p~K° and the last one is 
B~ — > ujK~ . We list in the following, the tree and penguin amplitudes which appear in 
the given transitions. 

For the decay B~ -> p°R- (a k = 32 in Eq. @), 
V2A r (a 1 ,a 2 ) = aifpF 1 (m 2 p )+a 2 fKAQ(rnj c ) , (63) 
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V2A*(a 3 , ■ ■ ■ ,a w ) = f p Fi(m 2 p )i^(a 7 + a 9 ) 

+ /xA (m^)|a 4 + aio - 2(a 6 + a 8 ) 

for the decay — >■ cufC - = 32 in Eq. (pip), 
\/2A^(ai,a 2 ) = a 1 f p F 1 (m 2 p ) + a 2 fKAo(m 2 K ) , 

V2A?(a 3 ,--- ,ai ) = / p Fi(mJ)|2(a 3 + a B ) + ^(a 7 + a 9 )| 
+ / i ^ (m^)<( -2(a 8 + a 6 ) 



A' 



(m u + m s )(m b + m u ) 



™A 



(m« + m s )(m b + m u ) _ 



+ (34 + aio 



for the decay B° -> p°if° (a* = 32 in Eq. (0)), 



V^AT(ai,a 2 ) = a x f p F x \ 



m 



P ) > 



>/2Af (a 3 , • ■ ■ ,a 10 ) = f p F x {m 2 p )^{a 7 + a 9 ) J 
+ /A-v4 (m|-)|a 4 - (2a 6 - a 8 ) 

for the decay B° -> cuiT (a fc = 32 in Eq. (0)), 
V / 2A^(ai,a 2 ) = a x f p F x (m 2 ) , 

\/2A^(a 3 , • ■ ■ ,a 10 ) = f p F x (rn 2 p )^2(a 3 + a 5 ) + ^(a 7 + a 9 ) 
+ /^^(m^)^ a 4 - (2a 6 



A 



(m s + m d )(m b + m d ) _ 



-aio 



m 



A 



_(ra s + m d )(m b + m d ) _ 



-aio 



for the decay B~ -> p^if (ar fc = 16 in Eq. (§!])), 



AT(ai,a 2 ) = a 2 f p F l (m 2 ) 



p) ' 



A pi. a 3, 



. aio) = fKA (m 2 K ) < a 4 - -a 10 - (2a 6 - 



m 



A 



(m s + m d )(m b + m d ) 
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for the decay B° -> p+K~ (a k = 16 in Eq. (§!])), 



A T Aa 1) a 2 ) = a 2 f K A (m 2 K ) 



(73) 



Ap(a 3 , ■ ■ ■ , a w ) = f K A Q (m 2 K ) < a 4 + a w - 2(a 6 + a 8 ) 



A' 



(m s + m u )(m b + m u ) _ 



• (74) 



Moreover, we can calculate the ratio between two branching ratios, in which the uncer- 
tainty caused by many systematic errors is removed. We define the ratio R as: 



R 



BR(5± -> p°^±) ' 
and, without taking into account the penguin contribution, one has, 



R 



BO 



1 + 



axfpF^m 



a2fKA (m 2 K ) 



1 + 



n 



pUJ 



s p - m l) + im u T u 



(75) 



-2 



(76) 



6.3 Numerical results 

The numerical values for the CKM matrix elements Vj ,F ', the p — to mixing amplitude 
Tipu,, and the particle masses mv,p, which appear in Eq. flBTJ), have been all reported in 
Section^. The Fermi constant is taken to be Gp = 1.166391x 10 -5 GeV~ 2 [|l^], and for the 
total width B decay, Vb(— 1/ t b) ; we use the world average B lifetime values (combined 
results from ALEPH, Collider Detector at Fermilab (CDF), DELPHI, L3, OPAL and 
SLAC Large Detector (SLD)) 



(77) 



t b o = 1.546 ± 0.021 ps 
r B+ = 1.647 ±0.021 ps 



To compare the theoretical results with experimental data, as well as to determine the 
constraints on the effective number of color, iV®, the form factors, and the CKM matrix 



parameters, we shall apply the experimental branching ratios collected at CLEO f32|| , 
BELLE []33], |34], |35| and BABAR 0, [I7| factories. All the experimental values are sum- 
marized in Table |7|. 

In order to determine the range of N^'* available for calculating the CP violating 
parameter, a, in B ±,Q — > p°K ±,(i , we have calculated the branching ratios for B^ 1 — > p°K ± , 
B ± -> p ± K°, B° -> p ± K^, B° -> p°K°, and B ± -> We show all the results in 

Figs. 0, [l^, [l^, and |18|, where branching ratios are plotted as a function of N%** for 
models (1) and (2) (different form factors are used in models (1) and (2)). By taking 
experimental data from CLEO, BABAR and BELLE Collaborations, listed in Table [7], 
and comparing theoretical predictions with experimental results, we expect to extract 
the allowed range of N^** in B — > pK and to make the dependence on the form factors 
explicit between the two classes of models: models (1, 3) and (5), and models (2) and (4). 
We shall mainly use the CLEO data, since the BABAR and BELLE data are (as yet) less 
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numerous and accurate. An exception will be made for the branching ratio — > u)K ± , 
where we shall take the BELLE data for our analysis since they are the most accurate 
and most recent measurements in that case. Nevertheless, we shall also apply all of them 
to check the agreement between all the branching ratio data. The CLEO, BABAR and 
BELLE Collaborations give almost the same experimental branching ratios for all the 
investigated decays except for the decay B~ — ► ujK~. In this later case, we observe a 
strong disagreement between all of them since they provide experimental data in a range 
from 0.1 x 10 -6 to 12.8 x 10 -6 . Finally, it is evident that numerical results are very sensitive 
to uncertainties coming from the experimental data and from the factorization approach 
applied to calculate hadronic matrix elements in the B — > K transition. Moreover, for 
B — > pK, the data are less numerous than for B — > pre, so we cannot expect to get a very 
accurate range of N*H. 

For the branching ratio B^ — > p°K ± (Fig. 14) we found a large range of values of 
Nl'* and CKM matrix elements over which the theoretical results are consistent with 
experimental data from CLEO, BABAR and BELLE. Each of the models, (1, 2, 3, 4) and 
(5), gives an allowed range of Aff ■ Even though strong differences appear between 
the two classes of models, because of the different used form factors, we are not able to 
draw strong conclusions about the dependence on the form factors. For the branching 
ratio B^ — > p ± K°, (Fig. f[5|) , BELLE gives only an upper branching ratio limit whereas 
BABAR and CLEO do not. Our predictions are still consistent with the experimental 
data for all models, for a large range of N%** . In this case, the numerical results for 
models (1) and (2) are very close to each other and, we need new data to constrain our 
calculations. 

If we consider our results for the branching ratio B° — > p^K T (plotted in Fig. [16]), 
there is agreement between the experimental results from CLEO and BELLE (no data 
from BABAR) and our theoretical predictions at very low values of N^** and the CKM 
matrix elements. All the models (1, 2, 3, 4) and (5), give branching values within the range 
of branching ratio measurements if iVJ is less than 0.07. The tiny difference observed 
between models (1) and (2) comes from the form factor A (k 2 ) (where A (k 2 ) refers to the 
B to p transition taken at k 2 = m 2 K ) since in that case, the amplitude computed involves 
only the form factor A (k 2 ). For the branching ratio B° — > p°K° shown in Fig. [17], neither 
CLEO, BABAR nor BELLE give experimental results. Nevertheless, from models (1) and 
(2), it appears that this branching ratio is very sensitive to the magnitude of the form 
factor Fi(k 2 ) (in our case, Fi{k 2 ) is uncertain because h\ = 0.360 or 0.762 in models (1) 
and (2), respectively) since the tree contribution is only proportional to F\. Moreover, 
from the range of allowed values of N%**, we can estimate the upper limit of this branching 
ratio to be of the order 20 x 10~ 6 . Finally, we focus on the branching ratio B — > uiK^ 
which is plotted in Fig. |18| for models (1) and (2). We find that both the experimental 
and theoretical results are in agreement for a large range of values of -AT® • But, the 
models (1) and (2) do not give similar results because the form factor F\, applied in these 
models, is very different in both cases. Moreover, the dependence of the branching ratio 
on the CKM parameters p and rj indicates that it would be possible to strongly constrain 
p and r\ with a very accurate experimental measurement for the decay B~ — > uoK~ . 

To remove systematic errors in branching ratios given by the B factories, we look 
at the ratio, R, between the two following branching ratios: BR(5° — > p^K^) and 



19 



BR(B ± -> p°K T ). The ratio is plotted in Fig. p function of N%" , for models (1) 
and (2) and for limiting values of the CKM matrix elements. These results indicate that 
the ratio is very sensitive to both N*" and to the magnitude of the form factors. The 
sensitivity increases with the value of N^f and gives a large difference between models 
(1,3) and (5) and models (2) and (4). We found that for a definite range of AT* , all 
models investigated give a ratio consistent with the experimental data from CLEO. It 
should be noted that R is not very sensitive to the CKM matrix elements. Indeed, if we 
only take into account the tree contributions, R is independent of the CKM parameters 
p and rj. The difference which appears comes from the penguin contribution and has to 
be taken into account in any approach since they are not negligible. 

We have summarized for each model, each branching ratio and each set of limiting 
values of CKM matrix elements, the allowed range of N%** within which the experimental 
data and numerical results are consistent. To determine the best range of Aff, we have 
to find some intersection of values of Nl** for each model and each set of CKM matrix 
elements, for which the theoretical and experimental results are consistent. Since the 
experimental results are not numerous and not as accurate as one would like, it is more 
reasonable to fix the upper and lower limits of AT^" which allow us the maximum of 
agreement between the theoretical and experimental approaches. By using the limiting 
values of the CKM matrix elements we show in Table |8|, the range of allowed values of 
Av with p—uj mixing. Even though in our previous study for B — > pn, we have restricted 
ourselves to models (2) and (4) rather than models (1, 3) and (5), here, we cannot exclude 
one of the models (1,2,3,4) and (5) due to the lack of accurate experimental data. We 
find that iV c e// should be in the following range: 0.66(0.61) < iV c e// < 2.84(2.82), where 
the values outside and inside brackets correspond to the choice k 2 jm\ = 0.3(0.5). Finally, 
if we take into account the allowed range of N^** determined for decays such as B — > pn 
and B — > pK we find a minimum global allowed range of N^** which should be in the 
range 1.17(1.12) < Ntf* < 1.63(1.77). 



7 Summary and discussion 

We have studied direct CP violation in decay process such as B ±,Q — > p°K ±, ° — ► 7i + ti~K ± ' 
with the inclusion of p — uj mixing. When the invariant mass of the 7r + 7r~ pair is in the 
vicinity of the uj resonance, it is found that the CP violating asymmetry, a, has a max- 
imum a max . We have also investigated the branching ratios B° — > p°K°, B° — ► p ± K Zf , 
B ± — ► p ± K°, B ± — > p°K ± , and B ± — > uK ± . From our theoretical results, we make com- 
parisons with experimental data from the CLEO, BABAR and BELLE Collaborations. 
We have applied five phenomenological models in order to show their dependence on form 
factors, CKM matrix elements and the effective parameter N^f in our approach. 

To calculate the CP violating asymmetry, a, and the branching ratios, we started 
from the weak Hamiltonian in which the OPE separates hard and soft physical regimes. 
We worked in the factorization approximation where the hadronic matrix elements are 
treated in some phenomenological quark models. The effective parameter, A^f , was used 
in order to take into account, as well as possible, the non-factorizable effects involved in 
B — > pK decays. Although one must have some doubts about factorization, it has been 
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pointed out that it may be quite reliable in energetic weak decays |J8 . 

With the present work, we have explicitly shown that the direct CP violating asymme- 
try is very sensitive to the CKM matrix elements, the magnitude of the form factors Ao(k 2 ) 
and Fi(k 2 ), and also to the effective parameter N^* (in order of increasing dependence). 
We have determined a range for the maximum asymmetry, a max , as a function of the 
parameter N^->, the limits of CKM matrix elements and the choice of k 2 jm\ = 0.3(0.5). 
For the decay B° — > -k + -k~K° and from all models investigated, we found that the 
largest CP violating asymmetry varies from +37%(+55%) to — 20%(— 24%). As regards 
B- -> 7T+7T-K-, one gets +49%(+46%) to -22%(-25%). For B ±fi -> tt+tt-R^ , the 
sign of a max stays positive as long as the value of N c is less than 2.7. In both decays, the 
ratio between asymmetry values which are taken at upper and lower limiting p and r\ val- 
ues is mainly governed by the term 1/ sin/3. It appears also that the direct CP violating 
asymmetry is very sensitive to the form factors at high values of N^'* . We underline that 
without the inclusion of p — uj mixing, we would not have a large CP violating asymmetry, 
a, since a is proportional to both sin 5 and r. We found a critical point for which sin<5 
reaches the value +1, but at the same time, r becomes very tiny. We emphazise that 
the advantage of p — uj mixing is the large strong phase difference which varies extremely 
rapidly near the uj resonance. In our calculations, we found that for B ±,Q — ► 7t + 7t~K ± >°, 
the sign of sin£ is positive until N*" reaches 2.69(2.65) when k 2 jm\ = 0.3(0.5). Then, 
by measuring a for values of N^'* lower than the limits given above, we can remove the 
phase uncertainty mod(7r) in the determination of the CKM angle 7. 

As regards theoretical results for the branching ratios B — > p Q K ± 1 B^ — ^ p ± K°, 
B° — > p ± K T , B° — > p°K° and B ± — > ujK ± , we made comparison with data from the 
CLEO (mainly), BABAR and BELLE (for B ± -> wK*) Collaborations. We found that it 
is possible to have agreement between the theoretical results and experimental branching 
ratio data for B ± -> p ^, B ± -> p ± K°, B ± -> wK*, B° -> p ± K^, and R. For 
B° — ► p°K°, the lack of results does not allow us to draw conclusions. Only an estimation 
for the upper limit (20 x 10~ 6 ) has been determined. Nevertheless, we have determined a 
range of value of A^f , 0.66(0.61) < N^** < 2.84(2.82), inside of which the experimental 
data and theoretical calculations are consistent. We have to keep in mind that, because of 
the difficulty in dealing with non-factorizable effects associated with final state interactions 
(FSI), which are more complex for decays involving an s quark, we have weakly constrained 
the range of value of N^** . 

^From the CP violating asymmetry and the branching ratios, we expect to determine 
the CKM matrix elements. In order to reach our aim, all uncertainties in our calculations 
have to be decreased: the transition form factors for B — ► p and B — > K have to be 
well determined and non-factorizable effects have to be treated in the future by using 
generalized QCD factorization. Moreover, we strongly need more numerous and accurate 
experimental data in B — » pK decays if we want to understand direct CP violation in B 
decays better. 
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ure captions 



Fig. |T] Tree diagram, for B decays. 

Fig. H QCD-penguin diagram, for B decays. 

Fig. |3] Electroweak-penguin diagram, for B decays. 

Fig. H Electroweak-penguin diagram (coupling between Z, 7 and W), for B decays. 

Fig. |5] B decays without (upper) and with (lower) p — u mixing. 

Fig. P CP violating asymmetry, a, for B° — > 7T + 7T~K , for k 2 /m\ = 0.3, for N^f = 
0.66, 2.69, 2.84 and for limiting values, max (min), of the CKM matrix elements for 
model (1): dot-dot-dashed line (dot-dash-dashed line) for Nl'* = 0.66. Solid line 
(dotted line) for N e J f = 2.69. Dashed line (dot-dashed line) for Nf f = 2.84. 

Fig. |^ CP violating asymmetry, a, for B° — > 7i + ti~K , for k 2 /m 2 = 0.5, for N^f = 
0.61, 2.65, 2.82 and for limiting values, max (min), of the CKM matrix elements for 
model (1): dot-dot-dashed line (dot-dash-dashed line) for N^f = 0.61. Solid line 
(dotted line) for iV c e// = 2.65. Dashed line (dot-dashed line) for iV c e// = 2.82. 

Fig. ^ CP violating asymmetry, a, for B~ — ► 7r + 7T~K~ , for k 2 /m 2 = 0.3, for N^f = 
0.66, 2.69, 2.84 and for limiting values, max (min), of the CKM matrix elements for 
model (1): dot-dot-dashed line (dot-dash-dashed line) for N^* = 0.66. Solid line 
(dotted line) for N e J f = 2.69. Dashed line (dot-dashed line) for N* ff = 2.84. 

Fig. P CP violating asymmetry, a, for B~ — > Tt + -n~ K~ , for k 2 /m 2 = 0.5, for N^f = 
0.61, 2.65, 2.82 and for limiting values, max (min), of the CKM matrix elements for 
model (1): dot-dot-dashed line (dot-dash-dashed line) for N^* = 0.61. Solid line 
(dotted line) for iV c e// = 2.65. Dashed line (dot-dashed line) for iV c e// = 2.82. 

Fig. H sin 5, as a function of N^f, for B° -> n+n-R , for k 2 /m 2 B = 0.3(0.5) 
and for model (1). The solid (dotted) line at sin 5 = +1 corresponds the case 
n pw = (—3500; —300), where p — uj mixing is included. The dot-dashed (dot-dot- 
dashed) line corresponds to il p£J = (0; 0), where p — w mixing is not included. 

Fig. |TT| sin5, as a function of N%ff, for B~ — > 7i + 7i~K~, for k 2 /m 2 B = 0.3(0.5) 
and for model (1). The solid (dotted) line at sin 5 = +1 corresponds the case 
n pw = (—3500; —300), where p — uj mixing is included. The dot-dashed (dot-dot- 
dashed) line corresponds to Ii puJ = (0; 0), where p — u mixing is not included. 

Fig. [12] The ratio of penguin to tree amplitudes, function of Nl*, for B° — > 

7T + 7T~K°, for k 2 /m 2 B = 0.3(0.5), for limiting values of the CKM matrix elements 
(p,fj) max(min), for Il puJ = (—3500; — 300)(0, 0), (i.e. with(without) p — uj mixing) 
and for model (1). Figure 12a (left): for Yl pu} = (0;0), solid line (dotted line) 
for k 2 /m 2 B = 0.3 and (p, 77) max(min). Dot-dashed line (dot-dot-dashed line) for 
k 2 /m 2 B = 0.5 and (p, 77) max(min). Figure 12b (right): same caption but for Yl puJ = 
(-3500; -300). 
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Fig. O The ratio of penguin to tree amplitudes, r, for B — > tt + tt K . Same 



caption for Figure 13a (left) and Figure 13b (right) as in Fig. [12 



Fig. |1|] Branching ratio for B — > p°K ± , for models 1(2), k 2 /m 2 B = 0.3 and limiting 
values of the CKM matrix elements. Solid line (dotted line) for model (1) and max 
(min) CKM matrix elements. Dot-dashed line (dot-dot-dashed line) for model (2) 
and max (min) CKM matrix elements. 



Fig. |15| Branching ratio for B — > p ± K°, for models 1(2), k 2 /m 2 B = 0.3 and limiting 
values of the CKM matrix elements. Solid line (dotted line) for model (1) and max 
(min) CKM matrix elements. Dot-dashed line (dot-dot-dashed line) for model (2) 
and max (min) CKM matrix elements. 

Fig. Branching ratio for B° —>■ p ± K T , for models 1(2), k 2 /m 2 B = 0.3 and limiting 



values of the CKM matrix elements. Solid line (dotted line) for model (1) and max 
(min) CKM matrix elements. Dot-dashed line (dot-dot-dashed line) for model (2) 
and max (min) CKM matrix elements. 



Fig. [17| Branching ratio for B° — > p°K°, for models 1(2), k 2 /m 2 B = 0.3 and limiting 



values of the CKM matrix elements. Solid line (dotted line) for model (1) and max 
(min) CKM matrix elements. Dot-dashed line (dot-dot-dashed line) for model (2) 
and max (min) CKM matrix elements. 

Fig. [18] Branching ratio for 5 ± — ► uK^, for models 1(2), k 2 /m 2 B = 0.3 and limiting 
values of the CKM matrix elements. Solid line (dotted line) for model (1) and max 
(min) CKM matrix elements. Dot-dashed line (dot-dot-dashed line) for model (2) 
and max (min) CKM matrix elements. 

Fig. 19 The ratio of the two pK branching ratios versus iVj" for models 1(2) and 
for limiting values of the CKM matrix elements: solid line (dotted line) for model 
(1) with max (min) CKM matrix elements. Dot-dashed line (dot-dot-dashed line) 
for model (2) with max (min) CKM matrix elements. 
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Table captions 



• Table [I] Wilson coefficients to the next-leading order. 

• Table ^| Effective Wilson coefficients related to the tree operators, electroweak and 
QCD-penguin operators. 

• Table |3] Values of the CKM unitarity triangle for limiting values of the CKM matrix 
elements. 

• Table |4] Form factor values for B — > p and B — > K at q 2 = 0. 

• Table [5] Maximum CP violating asymmetry a max (%) for B° — ► tt + 7t~K , for all 
models, limiting values of the CKM matrix elements (upper and lower limit), and 
for k 2 /m 2 h = 0.3(0.5). 

• Table Maximum CP violating asymmetry a max (%) for B~ — ► 7r + 7r _ i^ _ , for all 
models, limiting values of the CKM matrix elements (upper and lower limit), and 
for k 2 /m 2 b = 0.3(0.5). 

• Table [7] The measured branching ratios by CLEO, BABAR and BELLE factories 
for B decays into pK (f0~ 6 ). 

• Table || Best range of N^f determined for k 2 /m 2 = 0.3(0.5) and for B — > pK 
decays (upper). Also range of N^f determined previously for B — > p7r decays [09] 
(updated). Finally global range of N^f from both B decays (lower). 
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Figure 3: Electroweak-penguin diagram, for B decays. 




q 3 =q 2 



Figure 4: Electroweak-penguin diagram (coupling between Z, 7 and W), for B decays. 
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Figure 5: B decays without (upper) and with (lower) p — u> mixing. 
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Figure 6: CP violating asymmetry, a, for B° — ► 7r + 7r _ ^ , for k 2 /ml = 0.3, for = 
0.66, 2.69, 2.84 and for limiting values, max (min), of the CKM matrix elements for model 
(1): dot-dot-dashed line (dot-dash-dashed line) for N%f* = 0.66. Solid line (dotted line) 
for N e J f = 2.69. Dashed line (dot-dashed line) for N e J f = 2.84. 




Figure 7: CP violating asymmetry, a, for B° — > tt+tt - ^ , for k 2 /m 2 = 0.5, for = 
0.61, 2.65, 2.82 and for limiting values, max (min), of the CKM matrix elements for model 
(1): dot-dot-dashed line (dot-dash-dashed line) for N^f = 0.61. Solid line (dotted line) 
for N e J f = 2.65. Dashed line (dot-dashed line) for N e J f = 2.82. 
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Figure 8: CP violating asymmetry, a, for B~ — ► n + 7T~K~, for k 2 jm\ = 0.3, for iVJ = 
0.66, 2.69, 2.84 and for limiting values, max (min), of the CKM matrix elements for model 
(1): dot-dot-dashed line (dot-dash-dashed line) for N^* = 0.66. Solid line (dotted line) 
for N^f = 2.69. Dashed line (dot-dashed line) for iV c e// = 2.84. 
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Figure 9: CP violating asymmetry, a, for B 



7T + 7T K 



for k 2 /m 2 b = 0.5, for N e Jf 



0.61, 2.65, 2.82 and for limiting values, max (min), of the CKM matrix elements for model 
(1): dot-dot-dashed line (dot-dash-dashed line) for N^f = 0.61. Solid line (dotted line) 
for N e J f = 2.65. Dashed line (dot-dashed line) for Nf f = 2.82. 
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Figure 10: sin 5, as a function of N^f, for B° — > 7t + tt~K°, for k 2 /m 2 B = 0.3(0.5) and 
for model (1). The solid (dotted) line at sin 5 = +1 corresponds to the case Ii piJJ = 
(—3500; —300), where p — uj mixing is included. The dot-dashed (dot-dot-dashed) line 
corresponds to II pw = (0; 0), where p — uj mixing is not included. 




Figure 11: sin5, as a function of iV^ for B~ — > 7T + 7T~K~, for k 2 /m 2 B = 0.3(0.5) and 
for model (1). The solid (dotted) line at sin 5 = +1 corresponds to the case n pw = 
(—3500; —300), where p — to mixing is included. The dot-dashed (dot-dot-dashed) line 
corresponds to Ii puJ = (0; 0), where p — uj mixing is not included. 
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Figure 12: The ratio of penguin to tree amplitudes, T '. clS cL function of N%** , for S° — > 
7r^~7r~ 7 for k 2 /m 2 B = 0.3(0.5), for limiting values of the CKM matrix elements (p, rf) 
max(min), for Il puJ = (—3500; — 300) (0, 0), (i.e. with(without) p — u> mixing) and for model 
(1). Figure 12a (left): for fl puj = (0; 0), solid line (dotted line) for k 2 /m 2 B = 0.3 and (p, 77) 
max(min). Dot-dashed line (dot-dot-dashed line) for k 2 /m 2 B = 0.5 and (p, rf) max(min). 
Figure 12b (right): same caption but for Yl puJ = (—3500; —300). 
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Figure 13: The ratio of penguin to tree amplitudes, r, for B~ 
for Figure 13a (left) and Figure 13b (right) as in Fig. 12. 



7r + 7r K . Same caption 
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Figure 14: Branching ratio for B ± — > p°K ± , for models 1(2), k 2 /m 2 B = 0.3 and limiting 
values of the CKM matrix elements. Solid line (dotted line) for model (1) and max (min) 
CKM matrix elements. Dot-dashed line (dot-dot-dashed line) for model (2) and max 
(min) CKM matrix elements. Notation: horizontal dotted line: CLEO data; dashed line: 
BABAR data; dot-dashed line: BELLE data. 
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Figure 15: Branching ratio for B ± p ± K°, for models 1(2), k 2 /m 2 B = 0.3 and limiting 
values of the CKM matrix elements. Solid line (dotted line) for model (1) and max (min) 
CKM matrix elements. Dot-dashed line (dot-dot-dashed line) for model (2) and max 
(min) CKM matrix elements. Same notation as in Fig. [T4|, but only experimental upper 
limits are available. 
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Figure 16: Branching ratio for B° —>■ p ± i^ =F , for models 1(2), k 2 /m 2 B = 0.3 and limiting 
values of the CKM matrix elements. Solid line (dotted line) for model (1) and max (min) 
CKM matrix elements. Dot-dashed line (dot-dot-dashed line) for model (2) and max 
(min) CKM matrix elements. Same notation as in Fig. |14| . 
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Figure 17: Branching ratio for B° — > p°K°, for models 1(2), k 2 /m 2 B = 0.3 and limiting 
values of the CKM matrix elements. Solid line (dotted line) for model (1) and max (min) 
CKM matrix elements. Dot-dashed line (dot-dot-dashed line) for model (2) and max 
(min) CKM matrix elements. 
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Figure 18: Branching ratio for B ± — > uK^, for models 1(2), k 2 /m 2 B = 0.3 and limiting 
values of the CKM matrix elements. Solid line (dotted line) for model (1) and max (min) 
CKM matrix elements. Dot-dashed line (dot-dot-dashed line) for model (2) and max 
(min) CKM matrix elements. Same notation as in Fig. |14| . 




Figure 19: The ratio of two pK branching ratios versus N%** for models 1(2) and for 
limiting values of the CKM matrix elements: solid line (dotted line) for model (1) with 
max (min) CKM matrix elements. Dot-dashed line (dot-dot-dashed line) for model (2) 
with max (min) CKM matrix elements. Same notation as in Fig. [nj. 
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) for \x = 5 GeV 






Ci 


-0.3125 






c 2 


+ 1.1502 




c 3 


+0.0174 


c 5 


+0.0104 


c 4 


+0.0373 


c 6 


-0.0459 


c 7 


-1.050 x 10" 


b c 9 


-0.0101 


C 8 


+3.839 x 10" 


4 c w + 


1.959 x 10~ 3 



Table 1: Wilson coefficients to the next-leading order (see the reference in text). 





g2/ m 2 = 3 




q 2 /ml = 0.5 




C[ 


-0.3125 




-0.3125 




c> 2 


+ 1.1502 




+ 1.1502 






+2.433 x 10~ 2 + 1.543 x 10" 


-H 


+2.120 x 10~ 2 + 2.174 x 10" 


-H 




-5.808 x 10~ 2 -4.628 x 10" 


-H 


-4.869 x lO" 2 - 1.552 x 10" 


~H 




+1.733 x 10~ 2 + 1.543 x 10" 


-H 


+1.420 x 10~ 2 + 5.174 x 10" 


~H 




-6.668 x 10~ 2 -4.628 x 10" 


-H 


-5.729 x 10~ 2 - 1.552 x 10" 


-h 


c 7 


-1.435 x 10~ 4 - 2.963 x 10" 


-H 


-8.340 x 10~ 5 - 9.938 x 10" 


-H 


C's 


+3.839 x 10~ 4 




+3.839 x 10~ 4 




C 'f 


-1.023 x 10~ 2 - 2.963 x 10" 


-H 


-1.017 x 10~ 2 - 9.938 x 10" 


-H 


Cio 


+ 1.959 x 10" 3 




+ 1.959 x 10" 3 





Table 2: Effective Wilson coefficients related to the tree operators, electroweak and QCD 
penguin operators (see the reference in text). 





a 


P 


7 


i^Pmini Vmin) 


104°47 


19°32 


56°21 


(Pmini TJmax) 


93°13 


24°31 


62°56 


(yPmax j Tjmin) 


112°14 


21°20 


46°66 


(Pmaxi Tjmax) 


99°66 


26°56 


53°78 



Table 3: Values of the CKM unitarity triangle for limiting values of the CKM matrix 
elements. 





h A 


hi 


m Ao 


mi 


do(di) 


60(61) 


model (1) 


0.280 


0.360 


5.27 


5.41 






model (2) 


0.340 


0.762 


5.27 


5.41 






model (3) 


0.280 


0.360 


5.27 


5.41 






model (4) 


0.340 


0.762 


5.27 


5.41 






model (5) 


0.372 


0.341 






1.400(0.410) 


0.437(-0.361) 



Table 4: Form factor values for B — > p and B — > K at q 2 = (see the reference in text). 
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AgL = 0.66(0.61) Ngf = 2.84(2.82) 



model (1) 


Pmax ■> Vmax 


32(46) 


-14(-16) 


Pmin ■> Vmin 


25(33) 


-19(-22) 


model (2) 


Pmax ) 'Umax 


32(41) 


-6(-7) 


Pmin ■> Vmin 


27(30) 


-9(-10) 


model (3) 


Pmax i 'Umax 


32(45) 


-14(-16) 


Pmin -i Vmin 


25(33) 


-20(-23) 


model (4) 


Pmax ■> Vmax 


32(41) 


-6(-7) 


Pmin ■> Vmin 


27(30) 


-9(-10) 


model (5) 


Pmax ) Vmax 


37(55) 


-15(-17) 


Pmin ■> Vmin 


26(40) 


-19(-24) 



Table 5: Maximum CP violating asymmetry a max {%) for B° — > 7r + 7r if , for all models, 
limiting values (upper and lower) of the CKM matrix elements, and for k 2 jm\ = 0.3(0.5). 





KLL = 0.66(0.61) 


N±{ x = 2.84(2.82) 


model (1) 


Pmax -i Vmax 


47(45) 


-15(-17) 


Pmin -i Vmin 


34(35) 


-21(-23) 


model (2) 


Pmax j Vmax 


45(41) 


-11(-13) 


Pmin j Vmin 


33(32) 


-17(-18) 


model (3) 


Pmax j Vmax 


47(44) 


-15(-17) 


Pmin -i Vmin 


34(35) 


-20(-23) 


model (4) 


Pmax ) Vmax 


45(42) 


-12(-13) 


Pmin -i Vmin 


33(32) 


-17(-18) 


model (5) 


Pmax -i Vmax 


49(46) 


-17(-19) 


Pmin -i Vmin 


36(35) 


-22(-25) 



Table 6: Maximum CP violating asymmetry a max (%) for B — > 7r + 7r K , for all models, 
limiting values of the CKM matrix elements (upper and lower limit), and for k 2 /m 2 = 
0.3(0.5). 
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CLEO 


BABAR 


BELLE 




8.46+^ ± 1.8" (< 17) 1 


10 ± 6 ± 2* (< 29) 1 


< 13. 5 1 


p ± K° 






< 23.6 1 


p±KT 


± 2.8* (< 32) 1 




1 K 8+ 5 - 1 + 1 - 7 * 
id - -4.6 -3.0 


p°K° 








BR{p ± K^ 1 ) 
BR(p°K ± ) 


1.89 ± 1.41 








3.2±|;2±0.8" (< 7.9) 1 


1.4±J;g ± 0.3* 


9.2^j±1.0* 



Table 7: The measured branching ratios by CLEO, BABAR and BELLE factories for 
B decays into pK (10 -6 ) (see the reference in text). Experimental data*, fit* and 
upper limit*'. 



B -> pK 




W\ 




model (1) 




0.66;2.68(0.61;2. 


68) 


model (2) 




1.17;2.84(1.09;2. 


82) 


maximum ran} 


re 


0.66;2.84(0.61;2. 


82) 


minimum ranj 


je 


1.17;2.68(1.09;2. 


68) 




B — > pit 








model (2) 




1.09;1.63(1.12;1. 


77) 


model (4) 




1.10;1.68(1.11;1. 


80) 


maximum ran} 


?e 


1.09;1.68(1.11;1. 


80) 


minimum ran}; 


;e 


1.10;1.63(1.12;1. 


77) 




global range 




{ Are//} 




global maximum range 


0.66;2.84(0.61;2. 


82) 


global minimum range 


1.17;1.63(1.12;1. 


77) 



Table 8: Best range of Nl** determined for k 2 jm\ = 0.3(0.5) and for B — > pK decays 
(upper). Also range of N*ff determined previously for B — > pn decays |^| (updated). 
Finally global range of N^f from both B decays (lower). 
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